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Abstract 

By solution of the Bethe-Goldstone equation for the Cooper pairing problem, an 
approximate analytic relation is derived between coherence length £ and the binding 
energy of the Cooper pair. This relation is then qualitatively confirmed by numeri- 
cally solving the corresponding self-consistent gap equations, following the crossover 
from weak to strong coupling, in non-s-wave superconductors. The relation applies 
to non-conventional superconductors, and in particular to heavy Fermions and to 
high-T c cuprates. Utilizing in addition a phenomenological link between k-oT c and a 
characteristic energy e c = h 2 /2m* £ 2 , with m* the effective mass, major differences 
are exposed in the functional relation between ksT c and e c for s-wave materials 
and for non-conventional superconductors. The relation between critical tempera- 
ture and e c thereby proposed correctly reflects the qualitative properties of heavy 
Fermion superconductors. 
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1 Introduction 



In earlier work [1,2], we have discussed the possible correlation between the 
critical temperature T c in anisotropic superconductors and a 'natural' en- 
ergy scale e c ~ h 2 /m*l 2 involving the effective mass m* and some char- 
acteristic length l c . Here, by anisotropic superconductors we mean a super- 
conductor characterized by non-spherically symmetric pairing, giving rise to 
an anisotropic k-dependence of the gap energy in momentum space. Well- 
known instances of such materials are the heavy Fermion compounds, most of 
them being characterized by a p-wave order parameter [3,4], and the high-T c , 
cuprates, whose order parameter displays d-wave symmetry [5]. The most nat- 
ural physical choice for the characteristic length l c entering the definition of e c 
above was the coherence length £ [1]. In the case of anisotropic superconduct- 
ing materials in the presence of magnetic fluctuations, we later correlated the 
coherence length £ to the spin- fluctuation temperature T sf as k B T si ~ h 2 /m*^ 2 
[2]- 

The possible enhancement of the critical temperature T c in an anisotropic su- 
perconductor was already studied in the early work of Markowitz and Kadanoff 
[6] within a weak-coupling, BCS-like approach, and more recently revived in 
Refs. [7,8,9] as possibly relevant for the high-T c cuprates. A generalization 
of Markowitz and Kadanoff 's results to the strong and intermediate coupling 
regime however showed that in the strong coupling limit anisotropy is effec- 
tively averaged out, and the gap tends to become isotropic [10]. 

In the case of the high-T c cuprates, the strong coupling limit corresponds to 
the underdoped region of their phase diagram, which is usually interpreted in 
terms of a crossover between Bose-Einstein condensation of strongly coupled 
preformed pairs above T c , and weak coupling, BCS-like superconductivity in 
the overdoped regime [11]. On the other hand, the coherence length £ can 
be continuously connected to the characteristic size of the preformed bosonic 
pairs in the normal state of (underdoped) high-T c cuprates. Therefore, the 
coherence length may serve to parametrize the crossover from weak to strong 
coupling, with A; F £ decreasing in going from weak coupling, characterized by 
large superconducting pair fluctuations, to strong coupling, k-p being the Fermi 
momentum [12] (see also Refs. [13,14], and refs. therein). In the underdoped 
regime, a relatively short coherence length is consistent with the idea of pre- 
formed pairs localized in real space [15,16,17,18]. In such a limit, it is interest- 
ing to investigate how the internal structure of these preformed pairs is related 
to the overall symmetry of the many-body order parameter [19]. 

Here, we derive an explicit expression for the coherence length in terms of the 
pair binding energy for the Cooper problem in anisotropic superconductors. 
Such an expression [see Eq. (13) below] explicitly contains the quantum num- 
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ber i of the pair relative angular momentum, which is usually employed to 
parametrize the anisotropic character of the order parameter, I = 0, 1,2 cor- 
responding to s-, p-, and d-wave symmetry, respectively. While this expression 
correctly reduces to the standard one for isotropic s-wave superconductors, in 
the case £ > it agrees qualitatively with the phenomenological dependence of 
kftTc on the characteristic energy e c , proposed in Ref. [1] for the heavy Fermion 
compounds as well as for the high-T c cuprates. Such a dependence of the char- 
acteristic energy for superconductivity on e c is then qualitatively confirmed 
by numerically solving the self-consistent gap equations for the maximum gap 
at T = in the case of an anisotropic superconductor in the crossover be- 
tween the weak- and strong-coupling limits, as a function of the dimensionless 
crossover parameter now employing a more general definition of the co- 
herence length £ for anisotropic superconductors. 

The paper is organized as follows. In Sec. 2, after briefly reviewing the Bethe- 
Goldstone equation for the Cooper problem in the case of isotropic supercon- 
ductors, we generalize and solve it for non-s-wave superconductors, character- 
ized by a superconducting instability in the (£, m) channel of relative angular 
momentum quantum numbers. Our results are discussed in relation to the 
earlier phenomenological findings in Ref. [1]. In Sec. 3 we analyze the gap 
equations arising from the corresponding many-body problem at the mean- 
field level. Although we now have to resort to numerical integration (at least 
in the non-s-wave case) , our results qualitatively confirm the approximate re- 
lation between characteristic energies derived in Sec. 2, also in the crossover 
from weak to strong coupling. In Sec. 4 we eventually summarize and propose 
some directions for future work. 



2 Bethe-Goldstone equation for non-s-wave superconductors 

The Bethe-Goldstone equation for the Cooper problem [20] in momentum 
space reads: 

(^-2£ k )^k = E^k^k', (1) 

k' 

where ^k is the Fourier transform of the pair wave-function with wave-vector 
k (here, we assume a spin-singlet state, with zero total momentum), £k = 
h 2 k 2 /(2m*) — /i is the energy of a single electron with respect to the chemical 
potential /x, m* is the effective mass, e is the binding energy of the electron 
pair, and \4k' is the Fourier transform of the electron-electron interaction 
(see also Ref. [21] for a pedagogical review). In the weak-coupling limit, we 
may safely assume fi = Ef at T = 0, with ep the Fermi energy. However, 
we anticipate that this identification will be superseded in Sec. 3, where the 
crossover from weak to strong coupling will be addressed more in detail. In 
the case of anisotropic superconductors, we adopt the spirit of the Anderson- 
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Brinkman- Morel model for p-wave superfluidity in 3 He [22], and expand the 
electron-electron interaction in spherical harmonics around the Fermi surface 

as 

Vkk' = -7rE E VWMk)^ m (k'), (2) 

" e=om=-e 

for |£k|, l^k'l < A, and zero otherwise, where k is the unit vector pointing 
along the direction of k, f2 denotes the volume of the system, and A is an 
energy cut-off, characterized by the nature of the interaction. In the case of 
conventional superconductivity (s-wave pairing, or I = 0), it would be natural 
to identify such an energy scale with the Debye energy, as in BCS theory. 

The use of spherical harmonics to expand the electron-electron interaction in 
anisotropic pairing superconductors characterized by a spherical Fermi surface 
was earlier considered by Markowitz and Kadanoff [6] in the weak-coupling 
regime (see also Refs. [8,9]), and later by Combescot [10] in the strong and 
intermediate coupling regime. A related model of exotic Cooper pairing with 
finite angular momentum has been discussed also in Ref. [23], where rotational 
symmetry breaking (I > 0) is due to the interplay between the finite range of 
the attractive potential and the interelectronic average distance. Within such 
model, the authors of Ref. [23] also derive the ^-dependence of the critical 
temperature in the Bose- Einstein limit. In the more realistic case of a non- 
spherical Fermi surface, spherical harmonics are naturally replaced by Allen's 
Fermi surface harmonics [24], which have been used by Whitmore et al. [7] 
in extending the results of Ref. [6] within the framework of Eliashberg equa- 
tions for T c . Moreover, spherical harmonics afford a natural classification of 
anisotropic superconductors (such as heavy Fermion compounds [3,4] as well 
as high-T c cuprates [5]) in terms of their pairing symmetry. 

Within the weak-coupling approximation, the largest attractive coupling con- 
stant Vi m = V > in Eq. (2) gives rise to a pairing instability in the (£, m) 
channel. In the following, we shall neglect other instabilities, which may arise 
well below T c , corresponding to mixed symmetry pairing. This amounts to 
retaining the (£,m) term only in the expansion Eq. (2). In this case, Eq. (1) 
has a solution given by 

belonging to the eigenvalue 

£ "- 2Aexp ( - wy ■ w 

where N(0) is the density of states at the Fermi level, and a is a normalization 
constant. Such a solution corresponds to a bound state (e < 0), and a further 
mean-field analysis of the many-electron problem (see also Sec. 3) shows that it 
indeed corresponds to a superconducting state, characterized by a gap function 
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at T = 0: 

A t = 2Are- 1 /^°)^, m (k), (5) 
having the symmetry of the attractive channel under consideration, where: 

lnT = -J ^ k |F £m (k)| 2 ln|F £m (k)|, (6) 

the integration being carried over the unit sphere [25]. The anisotropic k- 
dependence of the pair wave-function ^ k in Eq. (3) in the case £ > provides 
interesting information on the internal structure of a Cooper pair and its 
connection with the overall symmetry of the many-body gap function. 

Given the pair wave-function the coherence length £ is naturally defined 
by 



p 2 = Sk |Vk^k| m 

Passing to the continuous limit, in the isotropic, s-wave case (£ = 0), one 
obtains [21] 

e - ^ (8) 

where x = |e|/(2£p) measures the binding energy of a pair in units of the 
energy of an unbound pair. Apart from a numerical factor, from Eq. (8) one 
thus recovers the correct order of magnitude relation between the critical tem- 
perature, the Fermi velocity Vp = h^d^/dk, and the coherence length: 

k B T c ~ |e| ~ (9) 
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In the anisotropic case, for a pairing instability in the (£, m) channel, with 
the pair wave- function given by Eq. (3), the denominator in Eq. (7) is easily 
integrated in the continuous limit as 



4tt Jo (e-2E) 

8716 ' 



(10) 



where N(E)/4tt is the density of states per unit solid angle, and use has been 
made of the normalization condition of the spherical harmonics. In the case 
£ 7^ 0, the anisotropic k dependence of the pair wave-function Eq. (3) gives 
rise to two contributions in the numerator of Eq. (7), according to the chain 
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rule 

Here, we have made use of the decomposition Vk = & _1 V k + \td/dk, where 
V k denotes the angular part of the gradient operator in momentum space, 
and of the fact that £ k depends only on k = |k|. Passing to the continuous 
limit, integrating separately over the angles as in Eq. (10), and making use of 
the identity 



J ^ k |v t F fet (k)| 2 = - 1 dQ k y; m (k)v?y, m (k) 

= t(t + i), (12) 

which follows from a variant of the Green's formula over the unit sphere, one 
eventually obtains 



h 2 



4 J_ 1(1+1) ( _ x\nx 
3 x 2 1 + x \ 1 + x 



(13) 



Equation (13) implicitly relates the binding energy |e| of a Cooper pair to its 
characteristic size £ for anisotropic pairing superconductors. 

In the case of non-s-wave superconductors, Equation (13) manifestly includes 
an 'orbital' contribution, proportional to £(£+1), as earlier surmised in Ref. [1] 
on the basis of phenomenological considerations. There, we proposed the exis- 
tence of a phenomenological relation linking /cb^c to the characteristic energy 

£ < = (14) 

in the case of anisotropic superconductors. That the effective mass m* should 
enter inversely in determining the scale of k^T c was earlier recognized by Ue- 
mura et al. [26], who did not, however, include the coherence length in their 
analysis [12]. By comparing the experimental values for several heavy Fermion 
compounds (p-wave superconductors, £ = 1) as well as high-T c superconduc- 
tors (d-wave superconductors, £ = 2), in the anisotropic case we found that 
ksT c = f(e c ) deviates from the square- root behavior, k-QT c oc ^ l /e^, that is 
easily derived from Eqs. (9) and (14) for isotropic, s-wave superconductors. 
In particular, in the heavy Fermion case, we noticed a large initial slope in 
f(s c ) for e c = 0, and a tendency of such function to approach saturation, as 



In view of the fact that k^T c ~ \e\ ~ xep [21], Eq. (13) implicitly defines 
the generalization to the anisotropic case of the relationship between k^T c 
and e c we were looking for. Equation (13) correctly reduces to Eq. (8) in 
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the isotropic case (s-wave superconductors, £ = 0). In the case of non-s- 
wave superconductors (I > 0), Eq. (13) indeed increases with increasing e c , 
starting with a logarithmically infinite slope at e c = 0, and tending to a 
saturation value as e c — > oo. Such a limit is e.g. approached for /c F £ <C 1. 
Setting a = A/[3£(£ + 1)] and performing an asymptotic analysis of Eq. (13), 
in the limit k F £ <C 1 (£ > 0) one obtains 

X= W[(l + a)/eY (15) 

where W(z) is Lambert's function [27], and in the limit of very large anisotropy 
(£»1): 

lim a; = W-^e -1 ) =3.591.... (16) 

Figure 1 shows the dependence of the dimensionless measure of the pair bind- 
ing energy x = \e\/2e F on e c /e F = (k F £)~ 2 , implicitly defined by Eq. (13) in 
the cases £ = 0, 1, 2, corresponding to s-, p-, and c?-wave superconductivity, 
respectively. Keeping fixed all other variables, one notices that for £ > £ , 
where 

h 2 3 
2m*$e F = AW 2 {e-i) = 9 ' 672 ' ' ' ' 
|e| increases as £ increases while, for £ < £ , \e\ decreases as £ increases, 
although without large deviations from the limiting value, Eq. (16). Therefore, 
at least for sufficiently weakly coupled superconductors (£ > £o), anisotropy 
enhances superconductivity, the degree of anisotropy being here parametrized 
by the order £ of the spherical harmonic modelling the k dependence of the 
order parameter. 



3 Mean-field analysis of the many-body problem and the crossover 
from weak to strong coupling 



The definition of the Cooper pair size £ we employ in Eq. (7) makes use of 
the pair wave-function ifj^ for the Cooper problem, in which the only many- 
body effect is that of forbidding the occupancy of states below the Fermi 
level, e F . Within BCS theory, i.e. at the mean field level, a self-consistent 
treatment of the superconducting instability affords a better definition of the 
pair wave-function, which in momentum space is given by [21] ip^ = Ak/ (2-Ek), 
where A k is the gap function, now depending on k as a vector, and E k = 
(££ + (Akl 2 ) 1 / 2 the upper branch of the quasiparticle dispersion relation. With 
this definition of ip^, when E k is allowed to vanish locally on the Fermi surface, 
as is the case for p- and rf-wave superconductors, it has been recently shown 
that nodal quasiparticles give rise to a logarithmically divergent contribution 
to £, as defined by Eq. (7) [28]. While such a drawback does not arise in 
s-wave superconductors [29], in the case of anisotropic superconductors the 
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coherence length must be defined in terms of the range in real space of the 
static correlation function for the modulus of the order parameter [30,28]. 
The fact that our result Eq. (13) provides a finite estimate for £ also in the 
anisotropic case clearly depends on our approximate choice for which solves 
only the two-body Cooper problem. 

In order to discuss the limits of validity of the approximate results derived in 
Sec. 2, we will then work out numerically the mean-field solution to the corre- 
sponding BCS Hamiltonian, both in the weak and in the strong-coupling limit. 
We start by reviewing the model and notations set out in Refs. [12,30,29,28] 
(see also Ref. [31]). 

We shall consider the following Hamiltonian for a superconducting system in 
three dimensions [32]: 

H = 6ccL C ka + ^kk'(i T cl k+tU C_ k / +tU C k / T , (18) 

k<7 kk'q 

where c kcr [c kcr ] is the creation [destruction] operator for an electron with wave- 
vector k and spin projection a G {Til} along a specified direction, £ k = 
h 2 k 2 /(2m*) — /j, is the dispersion relation for free electrons with effective mass 
m*, measured with respect to the chemical potential /i, and 

= -^ m (k)y, m (k') (19) 

is the projection of the electron-electron interaction along the (£, m) channel, 
as in Eq. (2), which we assume to be attractive (V > 0) [33]. Equation (19) 
generalizes to the non-s-wave case the contact potential in real space discussed 
e.g. by Marini et al. [29]. Standard diagonalization techniques then lead to 
the mean-field coupled equations for the gap energy A k = A Ye m (k) and the 
particle density n at T = (see Ref. [31] for the case 0): 



™ = §E^k, (20b) 

k 

where t> k = |(1 — £k/-Ek), together with u 2 , = 1 — t> k , are the usual coherence 
factors of BCS theory. 

Owing to our choice of a free particle dispersion relation (all band effects 
are embedded in a single parameter, namely the effective mass m*), and of a 
contact potential in real space, the sums over k in Eq. (20a) above give rise 
to an ultraviolet divergence, which requires a suitable regularization. In three 
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dimensions, it is customary to intoduce the scattering amplitude a s [30,29,34], 
which for our anisotropic interaction reads: 




s 



(21) 



Subtracting Eq. (21) from Eq. (20a) one has: 




s 



(22) 



Following Ref. [29], we render Eqs. (20) and Eq. (22) dimensionless, by intro- 
ducing the dimensionless quantities 



and the Fermi energy e F = h 2 kp/(2m*) = h 2 (3ir 2 n) 2 / 3 /(2m*). In particular, 
Marini et al. [29] observe that x = /i/A can be used as a parameter for 
the crossover between the strong-coupling, Bose-Einstein (BE) limit (x <C 0) 
and the weak-coupling, BCS limit (xq ^> 0). Moreover, on one hand, the 
expression resulting from the dimensionless version of Eq. (22) for (/cfO s ) _1 
can be inverted to obtain x as a function of the dimensionless scattering length 
kpa s . On the other hand, one may alternatively use kp^ as the independent 
variable in place of kpa s [29]. Indeed, it was earlier recognized by Pistolesi and 
Strinati [12] (following the seminal work of Nozieres and Schmitt-Rink [15]) 
that a natural variable which can be used to follow the crossover from strong- 
coupling to weak-coupling superconductivity is the product /cf£ of Fermi wave- 
vector kp times the coherence length £ for two-electron correlation. In the BE 
limit, electrons are expected to bind in quasi-bound pairs localized in real 
space (Schafroth pairs [35]), thus realizing the condition kpC, <C 1, while the 
BCS limit corresponds to loosely coupled pairs, with /cp£ ^> 1, localized in 
momentum space close to the Fermi energy. 

A zero-temperature calculation of the coherence length for a three-dimensional, 
s-wave superconductor along the crossover between the weak- and the strong- 
coupling limits has been performed both numerically [30] and analitically [29]. 
The case of a two-dimensional, <i-wave superconductor has been discussed in 
Ref. [28], where a dispersion relation typical of the cuprate superconductors 
has been explicitly considered. As anticipated above, in the non-s-wave case 
(£ 7^ 0), the standard definition, Eq. (7), of the coherence length leads to un- 
physical divergences. One may still conveniently define a coherence length as 
the range in real space of the static correlation function X A (r) for the modulus 



h 2 k 2 1 




x = 



2m* A ' 




•o 
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of the order parameter [30,28] as 



r i = - l im l0 g^. (23) 

The Fourier transform A^(q) of such a function has been derived for an s- 
and a <i-wave superconductor in Refs. [30] and [28], respectively. In the case 
of our anisotropic interaction, Eq. (19), it reads: 



1 

2fi 



Em 



i 

E 



|2 (' U k+q/2^k-q/2 ~ ^k+q/2^k-q/2) Z 
^k+q/2 + -E'k-q/2 

£ + £--|A + ||A_|\ 



E+ + E_ 



E+E— 



(24) 



where E = E k , Y = Y^ m (k), £±, A±, £± are calculated at momenta k ± q/2, 
respectively, and use has been made of the gap equation, Eq. (20a), in going 
into the second line. 



First of all, since the summand in Eq. (24) depends only on k — |k|, q — |q|, 
and on the relative angle between k and q, passing to the continuum limit and 
transforming the sums over wave- vectors into an integral, one has X^(q) = 
Xa(q)- Moreover, X&(q) is an even function of its argument. Back to real 
space, one analogously has X^(r) = Xa(t). Then, the asymptotic behaviour 
of Xa(t) at large distance r will be governed by the behaviour of its Fourier 
transform X&(q) at small wave- vector q. In particular, assuming the expansion 
[30,28]: 

X A 1 (q)=a + bq 2 + 0(q A ), (25) 
it is straightfoward to show [30] that 

e = -. (26) 

a 

Such a definition of the coherence length is now consistent also for nodal 
superconductors [28]. It reduces to Eq. (7) in the weak-coupling limit for the 
s-wave case, and applies also to the strong-coupling regime, provided b > 
(given that a > 0, identically, provided A ^ 0), i.e. provided that X^{q) has 
its absolute minimum at q = [36]. 

We have numerically solved the gap equations, Eqs. (20), and evaluated the 
dimensionless coherence length kp^ according to this more general definition, 
Eq. (26), for the anisotropic potential Eq. (19), with £ = 0, 1,2, correspond- 
ing to s-, p-, and <i-wave pairing, respectively. Here, a convenient measure 
of the characteristic energy for superconductivity may be taken as the maxi- 
mum gap A max , where A max oc k^T c holds also for anisotropic superconductors 
[10]. Such a quantity is readily extracted from the solution of the gap equa- 
tions, Eqs. (20). In Figure 2, we plot the dimensionless characteristic energy 
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2A max /e F versus e c /e-F = (&f£)~ 2 i n the cases i = 0,1,2. As k F t; decreases 
[(^fO~ 2 increases], one crosses over from the weak-coupling, BCS limit into 
the strong-coupling, BE limit [12]. Despite a mean- field solution of the self- 
consistent gap equations has been now taken into account in the calculation, 
and a more general and consistent definition of the coherence length has been 
employed, our numerical results are in good qualitative agreement with the 
approximate results of Sec. 2, with 2A max /£p increasing as a function of e c /e-f 
with a steeper slope, as t increases from t = (s-wave) to £ = 2 (d-wave), 
all curves tending to saturation, at least within the numerically accessible 
range of the crossover parameter fc F £. These results are in agreement with the 
phenomenological plots correlating characteristic energies for both the heavy 
Fermion and cuprate superconductors in Refs. [1,2]. 



4 Summary and future directions 

We have studied the Cooper problem for an anisotropic superconductor char- 
acterized by an electron-electron interaction expanded in terms of simple 
spherical harmonics over the Fermi sphere. In the weak-coupling limit of a 
superconducting instability in the (£, m) channel, we have derived an ana- 
lytical expression for the relation between the pair binding energy and the 
correlation length, for arbitrary relative angular momentum quantum num- 
ber I. While such an expression correctly reduces to the standard one in the 
s-wave case (£ = 0), in view of the fact that k^T c scales with such a pair 
binding energy, in the case of non-s-wave superconductors (£ > 0) our ex- 
pression agrees qualitatively with the phenomenological correlation between 
k&T c and the characteristic energy h 2 /(2m*^ 2 ) earlier found in Ref. [1]. These 
results have been confirmed by a numerical solution of the self-consistent gap 
equations in the crossover between the weak- and the strong-coupling limits, 
where a more general definition for the coherence length has been employed, 
applying to anisotropic superconductors. 

It may, in the future, prove of significance to include the effect of Coulomb 
interaction on the formation of Cooper pairs in superconducting assemblies 
along lines such as those laid down in Ref. [37] (see also Ref. [38]). The question 
arises then as to whether, in the future, it will be of significance in making the 
present study the basis of fully quantitative calculations, to generalize beyond 
the result here, and of course, of many earlier workers, that the Cooper pairs 
formed as a consequence of an attractive interaction correspond to a single 
value of the binding energy at a given temperature. There is some evidence 
(see, e.g., Refs. [39,40]) that some physical properties, and in particular specific 
heat and tunneling spectra of cuprate materials, may require generalization to 
Cooper pairs with a finite range of binding energies. We do not anticipate that 
the qualitative trends proposed in the present paper will be grossly affected. 
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Fig. 1. Dependence of the normalized Cooper pair binding energy |e|/2eF on the 
characteristic energy e c /ep = [h 2 /(2t?7*£ 2 )]/£f = (^f0~ 2 > as implicitly defined by 
Eq. (13), for £ = 0,1,2, corresponding to s-, p-, and d-wave superconductivity, 
respectively. See text for discussion. 
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Fig. 2. Dependence of twice the normalized maximum gap energy 2A max /eF on the 
characteristic energy £ c /£f = (&f£)~ 2 > f° r ^ = 0>1>2, corresponding to s-, p-, and 
d-wave superconductivity, respectively. As /cfC decreases [i.e., (&f£)~ 2 increases] one 
crosses over from the weak-coupling, BCS limit, to the strong-coupling, BE limit. 
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